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Abstract. We describe a quantum state tomography scheme which is applicable to a 
system described in a Hilbert space of arbitrary finite dimensionality and is constructed 
from sequences of two measurements. The scheme consists of measuring the various 
pairs of projectors onto two bases -which have no mutually orthogonal vectors-, the two 
members of each pair being measured in succession. We show that this scheme implies 
measuring the joint quasi-probability of any pair of non-degenerate observables having 
the two bases as their respective eigenbases. The model Hamiltonian underlying the 
scheme makes use of two meters initially prepared in an arbitrary given quantum state, 
following the ideas that were introduced by von Neumann in his theory of measurement. 
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1. Introduction 

The fundamental problem of inferring the initially unknown state of a quantum system 
from a set of measured quantities, i.e., the problem of quantum-state reconstruction, 
goes back to the early days of quantum mechanics, when it was known as Pauli problem 
PP. Aside from its fundamental importance, the ability to reconstruct and characterize 
quantum states has implications in various areas of technology and information sciences. 
Since a large amount of theoretical and experimental work has been devoted to this 
problem, we refer the reader to a representative number of contributions only and 
references therein [21E1H1E1E1IZ1IS1I2]. 

The main goal in formulating a tomographic scheme is to identify a set of 
measurable quantities that gives complete information about the state of the system. 
It was shown [3] that a complete set of mutually unbiased bases (MUB) is of this kind 
and could be used for state tomography with high efficiency However, the construction 
of a complete set of MUB is known for powers of primes dimensions only [10J. A 
possible alternative to MUB for state tomography is to use what are known as symmetric 
informationally complete probability-operator measurements (SIC POMs) [8J. These 
kind of measurements are known to exist in all dimensions d < 45 (with high numerical 
precision) [11]. So far, however, all experiments and even proposals for experiments 
implementing SIC POMs have been limited to the very basic quantum system, the two- 
level system (qubit) [12]. This is, in part, due to the fact that there is no systematic 
procedure for implementing SIC POMs in higher dimensions, in a simple experimental 
setup. 

In this contribution we identify sets of observables which, when measured in 
succession, provide complete information about the state of a quantum system described 
in a d- dimensional Hilbert space. The tomographic scheme consists of measuring the 
various pairs of projectors onto two bases which have no mutually orthogonal vectors, 
the two members of each pair being measured in succession. We formulate the scheme 
using the von Neumann (vNM) model for measurements extended to two successive 
measurements. We find that in the formulation the notion of finite-dimensional quasi- 
probability distribution (the analogue of the phase-space quasi-probability distribution) 
appears in a natural manner. 

The paper is organized as follows. In the next section we describe the dynamics of 
successive measurements of two arbitrary observables as an extension of the vNM model, 
assuming that the initial state of the two meters is described by an arbitrary density 
operator. In Sec. [3] we use this formalism to discuss the reconstruction scheme based 
on successive measurements. This procedure generalizes that of Ref. [2] to the more 
general case discussed in the present article. In Sec. H] we show that the notion of quasi- 
probability distribution and generalized transform of observables that was introduced 
in Ref. [H] can also be applied to the present general case. Finally, we present our 
conclusions in Sec. [5] 
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2. The successive measurements of two observables 

Consider the successive measurement of two observables: A is measured first and B 
later. For this purpose we generalize the standard vNM [151 IB], assuming two meters 
Mi (with canonical momentum and coordinate operators Pi, Qi, i = 1,2) which interact 
successively with the system of interest according to the interaction 

V(t) = e igi (t - h)AP 1 + e 2 g 2 (t - t 2 )BP 2 , (1) 

with < t\ < t 2 . The functions g\ and g 2 are normalized to 1 and have a (non- 
overlapping) compact support around the corresponding interaction times t\ and t 2 . 
We denote by (i = 1, 2) the strength of the interaction between the system and meter 
Mi. 

Before the first interaction, the system and the two meters are described by the 
density operators p s , pu^ and pm 2 , respectively, and the combined system by their tensor 
product p = PsPm 1 Pm 2 - After the second interaction, the combined system is described 

by m 

nn'mm' 

■ e- ieianA p Ml e ieia -' Pl ■ e- ie2bmp2 pM 2 e ie2b ™' p \ (2) 

where (B A) indicates that A has been measured first and B later. Of course, the 
whole process could be considered as one global measurement, leading to the density 
operator (j2J). Here, A and B are expressed in their spectral representation 

i = ^a n P a „, (3a) 

n 

B = Y,b m n m - (36) 

m 

The eigenprojectors F an and P& m correspond to the possibly degenerate eigenvalues a n 
of A and b m of B, respectively; they satisfy the orthogonality and completeness relations 

Pa„Fa„, =<WPa„, (4a) 
n 

and similarly for B. 

We now seek information on the system by observing the position-position 
correlation function of the two meters, (QiQ 2 ), where the average is taken over the 
state of Eq. (j2J). Under the assumption that the average initial position of the two 
meters vanishes, we obtain 

€\€ 2 

nm 

where 3? stands for the 'real part' and we have defined 

w£^ A \e 1 ) = £ A( £l (a n - O) tr(p s P an/ n m P Q J . (56) 
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The function A(ei(a n — a' n )) is explicitly given in Appendix A 

Consider again the same Hamiltonian of Eq. (pQ), but suppose that at the end of 
the measurement procedure (i.e., for t > t 2 ) we observe the momentum Pi of the first 
meter (instead of its position) and the position Q 2 of the second meter. The resulting 
momentum-position correlation function is 

— (Ag 2 ) { ^ A) = 2a 2 Pi ^ a nb m W^Z A \e 1 ) (6a) 

nm 

where 3" stands for 'imaginary part', ap ± is the second moment of Pi in the initial meter 
state, and we have defined 

w£t n A \ti) = E^( a « - OMp^nM- (66) 



The function A(ei(a n — a' n )) can be found in Appendix A as well. The quantity 
W^ B ^ A \e\) of Eq. fl5fe|) is in general different from W^ B ^~ (ei) of Eq. fl66|) . because 
A is in general different from A. For pure Gaussian states of the meters, A = A and 
therefore W = W. This was the starting point for the reconstruction scheme presented 
in [H]. In what follows we formulate a tomographic scheme in which this equality may 
not hold. 



3. State tomography scheme 

We now use the above formalism to describe a state tomography scheme. For this 
purpose we consider a d- dimensional Hilbert space and two orthonormal bases, whose 
vectors are denoted by \k) and respectively, with k, \i = l,...,d. Latin letters 
will be used to denote the first basis while Greek letters will be used for the second 
basis. We assume the two bases to be mutually non-orthogonal, i.e., (/c|/i) 7^ 0, Vfc,/x. 
This last condition implies that the two bases have no common eigenvectors and are 
said to be complementary [18]. The condition that the two bases have no common 
eigenvectors is equivalent to requiring that two observables having these bases as their 
eigenbases should never possess simultaneous definite values. An example of mutually 
non-orthogonal bases are two bases which are related by the Fourier transform. 
We now consider the following meters-system interaction 

V(t) = e l9l (t - t^A + e 2 g 2 (t - t 2 )¥^P 2 , (7) 

with < t\ < t 2 . Here, P^ = \k)(k\ and P M = are rank-one projectors onto the k- 

and /i-state of the first and second basis, respectively. The observable A appearing in 
the interaction of Eq. ([1]) is replaced here by the projector and B by the projector P^. 
Being projectors, these observables possess two eigenvalues: and 1. We denote by r 
and a the eigenvalues of P^ and P M , respectively, and the corresponding eigenprojectors 
by (P fc ) T and (P^. We have 

(Pfc)i = Pfc (8 a) 

(P fc ) = l-P fc = Pfc '' ( 86 ) 
fc'(^fc) 
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and similarly for (P M )i and (P jU )o- 

In the present case, Eq. fl5aj) for the meters position-position correlation function 
gives 

i 



1 



eie 2 



■(Q1Q2) 



T,a=0 



11 



(9 a) 



(96) 



In Eq. (J9a|), W ( 



ei) is the particular case of the quantity W^ A \e x ) of Eq. §R 



when A, £?, a n , and b m are replaced by P^, P M , r, and a, respectively, i.e., 

1 

W£m<-p*)( Ci ) = £ A(ex(r - r'))tr [p s (P,) r KP /i ) CT (P & ).] , 



r'=0 



and, in particular, 

W^"^ Pfc) (ei) = tr( Ps P fc P M P fe ) + A(ei) E tr(p s P fe ,P M P fe ) , 



(10) 



fc'(^fc) 



where we have used Eq. (jl.2d) of Appendix A Now we invert Eq. (ITTj) to obtain p s . 
This inversion, which was briefly indicated in Ref. [Hj, is valid under the more general 
situation contemplated here. To see this, we first write Eq. (11 ip as 



J2 G k'k(ei)(k\ Ps \k')(k'\iJ,) 



(n\k) 



where 



Gk>k( e i) — ^fe'fc + A(ei)(l ~ Sk'k)- 
From Eq. (|T2|) we obtain 

E /.....v W) = G fc / fc (e 1 )(fc|p a |fc') J 



so that 



and hence 



(Ail*) 



= E 



<^i£i).»l 



(12) 



(13) 



(14) 



(15a) 



P, = E W 



11 l e ly 



|fc> 1 



A(6l) 



(*l + l*> 



A(ei)<H*> 



(Ml 



(156) 

Eqs. fTT5fflj) and (pS5|> are the main result of this paper. They imply that V7c, /i, W / 1 ^ M ^ Pfc ^ 
of Eq. pip , just as p s , contains all the information about the state of the system. 
Therefore, inferring W 1 ^ >M<_Pfe ' ) Vfc, /i from the measurement outcomes is equivalent to 
the reconstruction of p s . Our aim is to show that the measured quantities, the position- 
position and momentum-position correlation functions, are informationally complete: 
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/p j p \ 

that is, one can reconstruct from these quantities. Note that neither the 

strength of the second interaction nor the state of the second meter enter Eq. ( 115&I) . 

/p j p \ 

We notice that the full complex quantity is needed for tomography. From 



the position-position correlation function of Eq. (I95j) we directly extract the real part 

/p ^ FY 1 ) 

of W£ M w . To find the imaginary part of Wn we measure the momentum-position 
correlation function, which in the present case is given by 

— (Ag 2 ) (P ^ Pfc) = 2a 2 P ^W^ ¥k \e 1 ), (16) 
eie 2 

where 

W^ ¥h \e 1 ) = tr(p s P fc P /t P fe ) + A(e x ) £ tr^P^P,) . (17) 

We note in passing that this equation can be inverted to write p s in terms of W, resulting 
in Eq. (11 5 aft with W and A are replaced by W and A. Though the function W is in 
general not equal to the function W, we prove in Appendix B that it contains all the 
information about the imaginary part of Wn , and therefore enables a complete 
state reconstruction. This completes our procedure. 

At first glance it seems that, in a d- dimensional Hilbert space, the present scheme 
for state reconstruction requires the measurement of the meters position-position and 
momentum-position correlations for the d 2 successive measurements of projectors, P& 
followed by P M , giving 2d 2 different measurements altogether. However, Hermiticity 
and the unit value of the trace of the density matrix p s impose d 2 + 1 restrictions 
among its matrix elements, so that p s can be expressed in terms of d 2 — 1 independent 
parameters. These restrictions eventually imply that only d 2 — 1 of these correlations are 
actually independent and thus the measurement of only d 2 — 1 correlations is required. 



Appendix C gives an application of the above formalism to the case of a two-dimensional 



Hilbert space, and shows explicitly how the matrix elements of p s can be expressed in 
terms of d 2 — 1 = 3 independent measurable correlations. 

Finally, we close this section with the following remarks. In the infinitely-strong 
coupling limit, t\ — > oo, A(ei) vanishes, and Wn contains information only about the 
diagonal elements of p s , as can be seen from Eq. ( TTTT) . In the other extreme of weak 
coupling, in particular in the limit when e\ — > 0, Wn contains the full information about 
the state of the system, cf. Eq. ( 1TTT) . This limit was the result presented in Ref. [13J. 
Therefore, to reconstruct a quantum state using the successive-measurement scheme it 
is better to perform a measurement with a weak coupling to the first meter rather than 
one with a strong coupling. 



4. A quasi-distribution and a generalized transform of observables 

From a conceptual point of view, one attractive feature of the present approach is related 
to the quantities Wn^ Fk \ei) that enter the state reconstruction formula, Eqs. 1 1 15 at) 
and (11551) . This quantity can be interpreted as a "joint quasi-probability distribution" 
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in the following sense. Let O be an observable associated with a (i-dimensional quantum 
system. Making use of Eq. (11 5 aft we can express its expectation value as 

tv(p s 6) = J2(k\Ps\k')(k'\6\k) =Y,W^ Fk \e 1 ) 0(a*,Mi), (18) 

kk' k[i 

where we have defined the "transform" of the operator O as 

0*1*0 (k'\d\k) 



0(ji, k,ei) = ^2 



(19) 

Eq. (118p has the structure of a number of transforms found in the literature, that express 
the quantum mechanical expectation value of an observable in terms of its transform and 
a quasi-probability distribution. For example, the Wigner transform of an observable 
and the Wigner function of a state are defined in the phase space (q,p) of the system, 
q and p labelling the states of the coordinate and momentum bases, respectively. In 
the present case, the transform ( TT9|) of the observable is defined for the pair of variables 
(//, k), /i and k labelling the states of each of the two bases. As Eq. (1T8|) shows, it is 
the quantity W^^^i^i) which plays the role of the quasi-probability for the system 
state p s , and is also defined for the pair of variables (p, k). It can be thought of as the 
joint quasi-probability of two non-degenerate observables, with the two bases being their 
respective eigenbases. Since any pair of mutually non-orthogonal bases can be used, we 
have a whole family of transforms that can be employed to retrieve the state. 

In the literature it has been discussed how Wigner's function can be considered as a 
representation of a quantum state (Ref. [17], Chs. 3 and 4), in the sense that i) it allows 
retrieving the density operator, and ii) any quantum-mechanical expectation value can 
be evaluated from it. Similarly, and for the same reasons, in the present context the 
quasi-probability W^' i *~ Fk \ei) can also be considered as a representation of a quantum 
state. 



5. Conclusions 



We discussed successive measurements as an alternative approach to realize 
informationally complete measurements on quantum systems. Here we considered a 
particular Hamiltonian model for successive measurements that involves the system 
proper and both meters, and is an extension to two meters of vNM of measurement. 

In the approach presented in this paper we considered, in the (i-dimensional Hilbert 
space of the system, two complete, orthonormal bases, assumed to be mutually non- 
orthogonal, also called complementary. The observables needed for the present scheme 
are the projectors onto the basis vectors of each one of these two bases. We then showed 
that the set of all pairs of successive measurements of such projectors, one for each basis, 
allows the complete retrieval of the system state (Eqs. (115 qp . (jl5fc|) ). We proved that 
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the scheme can be formulated for arbitrary states of the meters, and for an arbitrary 
strength of the meter-system interaction. 

We showed that this procedure can be interpreted as measuring the joint quasi- 
probability of pairs of non-commuting observables, in a way similar to the state 
reconstruction based on measuring the quasi-probability in phase space provided by 
the Wigner transform of the state. 

As a final note we wish to point out two possible extensions that might be of 
interest. One is the state tomography of a continuous- variable system using successive 
measurements. The second generalization concerns state tomography when both the 
system and the meters are described in a finite-dimensional Hilbert space. 
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Appendix A. The quantities X(/3) of Eq. (J56D and A(/3) of Eq. fl6fr|) 

We are using the notation (• • ■)m 1 = tr(pA/ 1 • • ■). Let us define 

g(P) = (e-^K, (1.1a) 
H(3) = ^(e-^Q 1 e-'^) Ml , (1.16) 



where f3 = ei(o n — a n >). Since (Oi)mi = 0, Eq. fll.lfcj) is well defined when (3 — > 0. The 
function A(/?) is defined as 

\(P) = g(P) + 2h(/3). (1.1c) 

The functions g(/3) and h(f3) satisfy the properties 

g*(P)=g(-0), 3(0) = 1, (1.2a) 

h*(0) = -h(-p), h(0) = - l -{Q x Pi + PiQi) Ml = 0, (1.26) 

and therefore 

A(0) = 1. (1.2c) 

In writing Eq. fll.2fej) we have assumed the natural condition that the current density at 
the point Q\ for the first meter prior to the measurement 

J(Qi) = ^-( P QiA + APqx)^ (1-3) 
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vanishes VQi (mi being the mass of the first meter, and Fq 1 = \Qi)(Qi\)- Thus 

J{Qi)dQi =—(P 1 ) Ml =Q, (1.4a) 
mi 

QiJ(Qi)dQ! = ^—{QiPi + PiQi) Ml = 0, (1.46) 

from which Eq. (11.26j) follows. This assumption does not affect the generality of our 
results and could be lifted in a straightforward manner. 

The function A(ex(a n — a' n )) is defined as follows. We first define 

m = y^-, (1.5«) 

where g(f3) is given in Eq. (11. lap . We are assuming that (A)a/i = 0, so that A(/3) is 
well defined when /3 — > 0: indeed, we find the series expansion 



M/?) = -< + z|(PiV + -- (i.o/o 



We then define A(/3) as 



A(« = || (1.6.) 



i-VlH^) «>+'"• (1.66) 



Appendix B. Construction of W\\ from WJV\\ and 
If we write 

Wf x ^ k) = x lik + iy llk , (2.1a) 
W^ h) = x ldt + iy lik , (2.1b) 

the correlation functions, Eqs. (l9fej) and (|T6|) . become 

(Qig 2 ) (p ^ Pfe) 



x Mfc , (2.2 a) 



(Aq 2 ) ( 



-Pit) 



2a^ fc . (2.26) 

Our aim is to express y^ in terms of the measured quantities. We go back to the 
expressions flEE]) and for ^["^(ei) and Wff^^fci). The quantities A(ei), 
A(ei) are known if the state of the measuring apparatus M\ is known [see Eqs. (11. lap . 
(fTT6|) . (fTTcl) ] and JTS^]; we write 

A(ei) = A r (ei)+iA < (e 1 ), (2.3a) 
A(e 1 ) = A r (e 1 )+iA i (e 1 ). (2.36) 
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On the other hand, the traces appearing in Eqs. (II ip and (j!7p are unknown; we write 
them as 

tr(p s P fc P M P fc ) = r% (2.4a) 

tr(p s F k ,F^ k ) = + is^. (2.45) 

Using Eq. (I15al) . r° fc of Eq. (I2.4al) can be written in terms of measured quantities only, 

as 

r°, k = \(k\p)\ 2 (k\p s \k) = \(k\^J2^- (2-5) 

A*' 

We introduce the definitions fiTTaL b). §L3a} b). <^Ma] b) and the result O in Eqs. (HI]) 
and (fT7|) and write the latter as 

a; Mfc = |(/c|/i)| 2 ^ x M / fc + A r r Mfc - AjS^, (2.6a) 

y/iit = Air^ + A r s Mfc , (2.66) 
= Air^ + A r s Mfc . (2.6c) 

For every pair of indices /i, k we now have a system of three linear equations in the 
three unknowns r M fc, s^k and y^k, which can thus be expressed in terms of the measured 
quantities y^ and the x mk of Eq. (\2.2al b). The result for y^ k is 

3{A(ei)A*(ei)} / , l2 v^ \ |A( ei )| 2 _ 

K{A(e 1 )A*(e 1 )} V ^ / K{A(e 1 )A*(e 1 )} 

We have thus achieved our goal of expressing W^ FtJ "^ Fk \ei), and hence p s of 
Eq. (115 M) . in terms of the measured correlations of Eqs. (l2.2al o). 



Appendix C. State reconstruction for a two-level system 

As an illustrative example of the general formulation of Sec. [3j we consider the state 
reconstruction of a two-level system. 

We take the projectors P& with k — 0, 1 and P M with p = ± as projectors onto 
the eigenstates of the Pauli operators a z and a x , respectively. The measured quantities 
are the correlation functions (QiQ2)^ Ftl '^ Fk ^ arid (PiQ2)^ F ^ Fk ^ ■ In the case of spin-1/2 
particles, the interaction of Eq. (JTj) could be realized by cascading two Stern-Gerlach 
apparatuses, equipped with additional constant electric fields, and one measures the 
particle position-position and momentum-position correlations in the z and x directions. 

For simplicity of the presentation of the example, we choose the case in which A(ei) 



and A(ei) defined in Appendix A are real. It can be shown that this can be achieved if 
the original state of the first meter is an arbitrary mixture of pure states which, in the 
coordinate representation, are real and of definite parity. As a result, from Eq. (12.71) we 
find that the are related to the measurable y^ k as 

y^k = ~~~ ~Ufj,k ■ (3.1) 

A(ei) 
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The (k\p s \k') = puk' elements of the density matrix, Eq. (115 aj) , can be expressed in terms 
of the 8 real quantities x^, V^k of Eq. (12. lap (or x^, Vfik, using (l3.ip ).The conditions 
of Hermiticity and unit trace of pkk 1 give 5 relations among the 8 correlations x^ and 
y^k, so that only 3 of them are independent. One can choose 

(0i4) (F+ ^ o) ,„„ v 

x+o = , (3.2a) 



2/-o = — 2 — , (3.2c) 



(AQ 2 ) (p -^ Po) 

20-^6^2 

as the 3 independent correlations, and one finds that the remaining 5 can be expressed 
in terms of them as 



1 

x+i — 2 ~ x -0i 


(3.3a) 


1 

X-i = 2 — x +0' 


(3.36) 


2/+o = -y-o, 


(3.3c) 


y+i = y-o, 


(3.3d) 


y-i = -y-o ■ 


(3.3e) 



Finally, the p^y matrix elements can be written in terms of measured quantities as 

Poo = x+o + x-a, (3.4a) 

p lx = l- x +0 - x_ , (3.46) 
1 2 

^( e i) A(ei) 
Pio = Poi • (3-4aO 
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